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Abstract 

Electroosmosis  is  widely  observed  in  different  micro-scale  flows.  Flow  modification  due  to  electroosmosis  is  important  for  applications  related 
to  fuel  cells  such  as  micro  direct  methanol  fuel  cells  and  micro-scale  devices.  Mechanism  for  delivering  methanol  into  the  micro-fuel  cell  is  very 
important  process  affecting  its  power  density  and  often  they  need  to  be  stand  alone  systems  with  no  external  pumps  or  other  ancillary  devices  to 
eliminate  the  parasitic  power  loss  accounting  from  fuel  feed.  In  this  study,  numerical  simulation  has  been  carried  out  for  the  electroosmotic  effect  on 
pressure-driven  flows  in  a  serpentine  microchannel  whose  side  walls  are  subjected  to  variable  zeta  potentials.  It  is  observed  that  for  a  non-uniform 
zeta  potential,  the  secondary  vortex  pair  and  their  strength  of  motion  generated  by  electroosmosis  is  modified  asymmetrically,  corresponding  to 
the  degree  of  the  variation  in  and  the  magnitude  of  the  zeta  potential.  It  is  also  found  that  the  flow  profiles  are  quite  different  for  variable  zeta 
potentials  as  compared  to  a  constant  zeta  potential  applied  on  the  channel  walls.  The  flow  phenomena  at  the  bend  of  the  serpentine  channel  has 
also  been  investigated,  and  flow  control  is  possible  by  regulating  the  applied  zeta  potential  on  the  channel  walls.  It  is  observed  that  vortices  shift 
in  a  clockwise  direction  as  the  applied  zeta  potential  is  increased  on  the  walls  of  the  serpentine  channel.  Formation  of  additional  vortices  has  also 
been  observed  for  very  low  values  of  applied  zeta  potential. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  convective  transport  in  micro  domain  is  of  great  impor¬ 
tance  as  it  applies  to  various  application  such  as  Micro  Electro 
Mechanical  Systems  (MEMS)  devices,  micro-heat  exchangers, 
micro-fuel  cells,  etc.  [1].  Electroosmosis  is  the  bulk  movement 
of  liquid  near  a  stationary  surface  due  to  an  externally  applied 
electric  field  [2].  In  fuel  cells,  the  typical  fuel  such  as  hydrogen 
for  Proton  Exchange  Membrane  (PEM)  or  methanol  for  Direct 
Methanol  Fuel  Cell  (DMFC),  flows  through  the  serpentine  flow 
field  present  in  the  bi-polar  plates  of  the  cell.  Due  to  the  electric 
potential  generated  in  the  micro-fuel  cell,  a  potential  gradient 
is  developed  in  the  bi-polar  plate,  the  direction  of  which  is  per¬ 
pendicular  to  the  flow  direction.  The  presence  of  an  aqueous 
solution  in  the  flow  field  such  as  methanol  in  case  of  DMFC, 
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results  in  electroosmosis  driven  flow,  which  is  of  interest  for 
micro-fuel  cell  applications. 

There  have  been  a  few  studies  reported  in  the  field  of  elec¬ 
troosmosis  driven  flow  in  micro  DMFC.  Karimi  and  Li  [3]  have 
modeled  the  electroosmotic  flow  in  the  fuel  cell  membrane  in¬ 
cluding  the  electrokinetic  effect.  They  have  modeled  the  mem¬ 
brane  pore  to  determine  the  electroosmotic  flow  through  the 
membrane  for  different  geometrical  and  operating  conditions. 
They  have  reported  a  nonlinear  increase  in  electroosmotic  flow 
with  increase  in  pore  size  and  found  that  the  electroosmotic  drag 
coefficient  increases  linearly  with  pore  size. 

Several  researchers  have  investigated  the  effect  of  varying 
zeta  potential  on  electroosmotic  flow  (EOF)  in  microchannels 
[4-8].  Such  EOF  are  in  the  direction  of  the  pressure  driven 
flow.  The  characteristics  of  EOF  in  a  microchannel  depend 
upon  the  nature  of  the  zeta  potential,  i.e.,  whether  it  is  uniform 
or  nonuniform.  Fu  et  al.  [4]  used  Nemst-Planck  equation  along 
with  a  Navier-Stokes  solution  to  model  the  EOF  that  occurs 
when  a  step  change  in  zeta  potential  is  applied.  The  results 
indicate  that  a  step  change  in  zeta  potential  causes  a  significant 
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Nomenclature 

Co  ionic  concentration  in  bulk  solution 
Dh  hydraulic  diameter  of  channel  cross  section 
E  electric  field  intensity 

F  Faraday’s  constant 

Hc  channel  depth 

L  half  length  of  section  of  straight  segment  of  chan¬ 

nel 

Rc  radius  of  curvature  of  channel  axis  at  U-bend 

Re  Reynolds  number 

Ru  universal  gas  constant 

u  velocity  component  in  v-direction 

v  velocity  component  in  y-direction 

Viniet  velocity  of  fluid  at  inlet 

w  velocity  component  in  z-direction 

Wc  channel  width 

z<?  valence  number  of  univalent  fluid 

z+/-  valence  number  of  the  positive/negative  ions  in 
the  fluid 

Greek  letters 

€  permittivity  of  fluid 

k  Debye  length 

/ 1  viscosity  of  fluid 

0  applied  potential 

0  potential  due  to  electrical  double  layer 
p  density  of  fluid 

Pe  charge  density 

£w  zeta  potential  or  wall  potential 


variation  in  the  velocity  profile  and  also  in  the  pressure 
distribution. 

Micro-scale  debris  and  manufacturing  irregularities  during 
fabrication  of  microchannels,  and  adsorption  of  organics  during 
analysis  of  microchannels  produce  nonuniform  zeta  potential 
distributions  along  the  channel  walls.  Herr  et  al.  [5]  have  shown 
the  influence  of  zeta  potential  on  the  velocity  profile  and  sample 
dispersion  rate  for  electroosmotic  flow  in  cylindrical  capillar¬ 
ies  with  nonuniform  zeta  potential  distribution.  Similarly,  Chen 
et  al.  [6]  investigated  numerically  the  effects  of  step  change  in 
zeta  potential  in  a  cylindrical  microchannel  on  electroosmotic 
flow.  It  is  observed  that  the  variable  zeta  potential  along  the  flow 
direction  generated  regions  of  positive  and  negative  pressure 
gradients  in  the  flow  field.  Lee  et  al.  [7]  studied  electroosmotic 
flow  in  a  cylindrical  microchannel  with  non-uniform  zeta  poten¬ 
tial  distribution.  They  showed  formation  of  distorted  electroos¬ 
motic  velocity  profiles  and  flow  circulation  resulting  from  the 
axial  variation  of  the  zeta  potential.  However,  these  studies  are 
related  to  straight  microchannels  only. 

Zhang  et  al.  [8]  studied  the  two-dimensional  flow  pattern 
in  microchannels  with  large  aspect  ratios  with  heterogeneous 
zeta  potentials,  and  discuss  the  formation  of  secondary  EOF’s 
generated  by  wavelike  zeta  potential.  Their  results  show  that 


the  heterogeneous  zeta  potentials  could  generate  complex  flow 
patterns  and  enhance  mixing.  However,  their  study  may  not  be 
representative  for  microchannels  with  nearly  square  cross  sec¬ 
tion. 

Souders  et  al.  [9]  have  studied  combined  electroosmotic  and 
pressure  driven  flows  in  a  three-dimensional  microchannel  with 
walls  covered  by  charged  bands.  Under  the  influence  of  elec¬ 
troosmotic  flow  only,  it  is  found  that  there  is  a  significant  amount 
of  fluid  folding  and  stretching.  By  adding  pressure  gradient 
against  electroosmotic  flow,  they  showed  an  increase  in  the  fluid 
particle  residence  time.  The  helical  flow  structure  is  also  ob¬ 
served  for  this  combined  electroosmotic  and  pressure  driven 
flow. 

Fuel  feed  into  the  micro-fuel  cell  is  an  important  process 
affecting  its  power  density  and  often  there  is  a  need  to  mini¬ 
mize  power  loss  accounting  from  fuel  delivery.  Buie  et  al.  [10] 
presented  a  novel  design  and  performance  of  a  planar  silicon 
electroosmotic  pump  (EO)  for  methanol  fuel  delivery  in  DMFC. 
Though  the  EO  pump  could  not  deliver  enough  pressure  head  for 
fuel  delivery,  nevertheless  it  was  interesting  to  note  that  the  mi¬ 
cro  DMFC  has  shown  to  have  improved  polarization  and  power 
density  characteristics  with  the  EO  pump.  The  study  in  particu¬ 
lar  emphasizes  the  performance  enhancement  due  the  combined 
mechanism  of  electroosmotic  phenomena  and  diffusion  for  fuel 
delivery. 

More  recently,  Rawool  and  Mitra  [12]  studied  the  transport 
of  methanol  in  micro-fuel  cell  under  the  applied  potential  which 
is  perpendicular  to  the  direction  of  the  flow.  Such  conditions  are 
unique  to  fuel  cell  applications,  and  are  not  studied  prior  to  this 
work.  However,  in  their  work,  the  side  walls  are  subjected  to  a 
constant  zeta  potential,  which  is  not  always  the  case  in  micro¬ 
fuel  cell.  As  observed  in  the  existing  literature  [5,6],  the  zeta 
potential  tends  to  change  along  the  flow  field  as  the  fuel  gets  con¬ 
sumed  in  the  chemical  reaction  occurring  at  the  catalyst  layer, 
adjoining  the  flow  channels.  The  heterogeneous  zeta  potentials 
in  microchannels  can  also  be  achieved  through  chemical  modi¬ 
fications,  coupled  capillaries  and  integrated  systems  [8].  Hence, 
it  is  important  to  understand  the  flow  behavior  in  the  micro- 
channel  of  the  fuel  cell  under  variable  zeta  potential  conditions. 
To  this  end,  the  Poisson-Boltzmann  and  Navier-Stokes  equa¬ 
tions  are  solved  numerically  to  investigate  the  electroosmotic 
driven  flow  phenomena  for  a  range  of  operating  conditions. 

2.  Physical  problem 

In  Fig.  1  schematic  of  the  serpentine  microchannel  model 
is  shown.  The  channel  considered  here  is  a  three  dimensional 
serpentine  channel.  The  cross  sectional  area  of  the  channel  is 
100  [x  m  x  100  pim.  The  model  consists  of  a  single  ‘S’  shaped 
segment  of  the  serpentine  channel.  It  has  three  straight  lengths 
connected  by  two  semicircular  arcs.  In  Fig.  1  length  L  =  300  p,m 
and  the  radius  of  curvature  at  bend  is  Rc  =  200  p,m. 

3.  Governing  equations 

Electroosmosis  is  the  motion  of  fluid  caused  by  externally 
applied  electric  potential.  In  practice,  most  surfaces  contain  a 
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residual  negative  charge.  When  such  charged  surface  comes  in 
contact  with  ionizing  fluid,  the  ions  having  polarity  opposite 
to  the  surface  charge  are  attracted  towards  the  surface.  These 
ions  accumulate  near  the  surface,  forming  an  electric  double 
layer  (EDL).  In  this  layer  ion  concentration  changes  from  a 
maximum  value  (zeta  potential,  £w)  near  the  wall  to  a  neutral 
state  in  the  fluid  core.  The  thickness  of  this  layer  is  character¬ 
ized  by  Debye  length,  k,  which  is  defined  as  the  distance  from 
the  solid  surface  to  which  the  charge  drops  to  e~l  (37%)  of  its 
maximum  value.  The  typical  range  of  zeta  potential  and  Debye 
length  is  20-150 mV  and  10-9  to  10-6  m,  respectively  [11,13]. 
When  electric  potential  is  applied  across  a  capillary,  the  ions  in 
the  EDL  move  under  the  influence  of  the  electric  field.  This  ion 
movement  causes  bulk  fluid  motion  due  to  the  momentum  trans¬ 
fer  from  EDL  to  the  fluid  core.  Thus,  the  force  acting  on  ions 
is  accommodated  as  a  body  force  in  Navier-Stokes  equation, 
which  reads  as 

p(V  •  V)V  =  -VP  +  fiV2V  +  pEE  (1) 

where  pe  is  the  charge  density  and  E  is  the  electric  field.  For 
steady  flow  with  constant  properties,  the  continuity  equation  can 
be  written  as 

(2) 


The  electric  field  E  in  Eq.  (1)  has  two  contributions.  One  is  the 
field  due  to  applied  potential  (0)  for  which  governing  equation 
is 

V2</>  =  0  (3) 

The  second  contribution  is  the  field  due  to  EDL  which  is  gov¬ 
erned  by  the  equation 

vV=-—  (4) 

€ 

The  charge  density  is  given  by 

pE  =  FCoz-e(z~F*/R“T)  -  FC0z+  e (*+**/*« r)  (5) 

which,  for  a  symmetric  electrolyte  of  valence  number  ze  be¬ 
comes 

pE  =  2FC0  sinh  (6) 

where  0  is  the  potential  due  to  electrical  double  layer,  F  the 
Faraday’s  constant,  Ru  is  universal  gas  constant. 

3.1.  Boundary  conditions 

For  fluid  flow,  a  constant  velocity  corresponding  to  a  given 
Reynolds  number  is  specified  at  the  channel  inlet.  At  outlet,  a 


V-  V  =  0 
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constant  pressure  (atmospheric)  is  specified.  No  slip  boundary  condition  is  imposed  on  all  the  remaining  walls.  Following  is  the 
mathematical  representation  of  the  above  conditions: 


u  =  0 

v  =  constant 


for 


0  <  x  <  100  jut  m 

y  =  0 


w  =  0 


P  =  0,  for 


0  <  z  <  100  [x  m 
ARC  [x  m  <  x  <  ARC  +  100  jx  m 
y  =  0  p,  m 
0  <  z  <  100  jx  m 


a  —  0 
v  =  0  >  , 

w  =  0 


for  all  other  boundaries 


For  electric  field,  a  fixed  potential  is  specified  at  top  and  bottom  walls  of  the  channel,  as  shown  in  Fig.  1,  which  can  be  expressed 
as, 


0  =  0  V,  for 


z  =  0 

for  all  x  and  y 


0  =  O.1V,  for 


z  =  100  p  m 
for  all  x  and  y 


A  zeta  potential  and  Debye  length  are  specified  on  the  side  walls  of  the  channel  as  shown  in  Fig.  1 .  The  Debye  length  is  kept  constant 
at  1  x  10-7  m.  The  zeta  potential  is  varied  using  stepwise,  linear  and  constant  profiles,  to  investigate  the  effect  of  zeta  potential 
profile  on  the  electroosmotic  flow.  Also,  both  walls  are  subjected  to  different  zeta  potential  conditions.  In  total,  seven  different 
combinations  are  studied  and  it  is  expected  that  such  variations  in  zeta  potential  are  present  in  the  actual  operating  micro-fuel  cells. 
The  variable  zeta  potential  conditions  can  be  expressed  mathematically  in  the  following  manner: 


Case  1: 

Case  2: 


Constant  zeta  potential  of  —50  mV  on  both  the  side  walls. 

£w  =  —50  mV,  for  0  <  z  <  100  jx  m }  On  the  left  and  right  wall 

Stepwise  zeta  potential  on  right  wall  and  constant  zeta  potential  of  —50  mV  on  the  left  wall. 
t;w  =  —10  mV,  for  0  <  z  <  30  p,  m 

=  —50  mV,  for  30  p,  m  <  z  <  70  jx  m 

£w  =  —90  mV,  for  70  |x  m  <  z  <  100  |x  m  ^ 


>  On  the  right  wall 


£w  =  —50  mV,  for  0  <  z  <  100  jx  m}  On  the  left  wall 
Case  3:  Stepwise  zeta  potential  on  left  wall  and  constant  zeta  potential  of  —50  mV  on  the  right  wall. 
=  —50  mV,  for  0  <  z  <  100  jx  m}  On  the  right  wall 


>  On  the  left  wall 


t;w  =  —10  mV,  for  0  <  z  <  30  jx  m 
t;w  =  —50  mV,  for  30  p,  m  <  z  <  70  p,  m 

£w  =  —90  mV,  for  70  p,  m  <  z  <  100  |x  m 

Case  4:  Stepwise  zeta  potential  on  left  and  right  walls. 
t;w  =  —10  mV,  for  0  <  z  <  30  [x  m 
t;w  =  —50  mV,  for  30  |x  m  <  z  <  70  (x  m 

£w  =  —90  mV,  for  70  |x  m  <  z  <  100  \x  m  ^ 

Case  5:  Linear  zeta  potential  on  right  wall  and  constant  zeta  potential  of  —50  mV  on  the  left  wall. 


>  On  the  left  and  right  walls 


=  —0.01  —  0ffx  Z  V,  for  0  <  z  <  100  p,  m  |  On  the  right  wall 
=  —50  mV,  for  0  <  z  <  100  p,  m}  On  the  left  wall 
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Table  1 

Boundary  conditions  for  Cases  1-7 


Case 

Boundary  conditions 

Case  1 

Constant  zeta  potential  on  both  walls 

Case  2 

Stepwise  zeta  potential  on  right  wall  and  constant 
zeta  potential  on  left  wall 

Case  3 

Stepwise  zeta  potential  on  left  wall  and  constant 
zeta  potential  on  right  wall 

Case  4 

Stepwise  zeta  potential  on  both  walls 

Case  5 

Linear  zeta  potential  on  right  wall  and  constant  zeta 
potential  on  left  wall 

Case  6 

Linear  zeta  potential  on  left  wall  and  constant  zeta 
potential  on  right  wall 

Case  7 

Linear  zeta  potential  on  both  walls 

Table  2 

Properties  of  methanol 


Physical  property 

Value 

Dielectric  property 

Value 

Density  (kg  m-3) 

785 

Electric  conductivity 

4.4  x  lO"4 

(flm-1) 

Viscosity  (Pa  s) 

5.6  x  1CT4 

Relative  permittivity 

33.62 

of  inlet  velocity  corresponding  to  a  range  of  Reynolds  number. 
A  Reynolds  number  of  0.001  is  used  in  the  simulations.  Here, 
Reynolds  number  is  defined  as 

„  pVinlet  Dh 

li 


Case  6:  Linear  zeta  potential  on  left  wall  and  constant  zeta  potential  of  —50  mV  on  the  right  wall. 
Zw  =  —0.01  —  Z  V,  for  0  <  z  <  100  p,  m  j  On  the  left  wall 

t;w  =  —50  mV,  for  0  <  z  <  100  p,  m}  On  the  right  wall 


Case  7 :  Linear  zeta  potential  on  left  and  right  walls. 

fu,  =  -0.01  -  Xi°o-‘Z  for  o  <  Z  <  100  |X  m} 

These  boundary  conditions  are  also  summarized  in  Table  1  for 
easy  reference. 

4.  Numerical  technique 

The  channel  geometry  is  created  and  meshed  using  CFD- 
GEOM  modeler.  Structured  grid  is  used  for  simulations.  A  grid 
independence  study  is  carried  out  by  successively  refining  the 
grid,  increasing  the  number  of  elements  in  the  domain.  The  grid 
is  made  extra  fine  in  the  near  wall  region  where  maximum  gra¬ 
dients  occur.  For  a  grid  having  260,431  elements,  the  solution  is 
found  to  be  grid  independent.  Hence  this  mesh  size  was  used  in 
all  the  cases.  CFD  software  CFD-ACE+  is  used  for  the  simula¬ 
tions.  Upwind  scheme  is  used  in  the  CGS  +  Pre  solver  [14]  for 
velocity  and  electric  field  while  AMG  solver  is  used  for  pres¬ 
sure  correction.  Parametric  solver  is  used  to  solve  for  values 


On  the  left  and  right  wall 

where  Dh  is  the  hydraulic  diameter  of  the  channel,  defined  as 
D  =  4 (Wc  x  Hc) 
h  2(WC  +  Hc) 

Here  Wc  is  the  width  of  the  channel  and  Hc  is  the  depth  of  the 
channel  as  shown  in  Fig.  1 .  The  effect  of  zeta  potential  profile  on 
the  velocity  profile  in  the  channel  is  studied.  The  physical  and 
dielectric  properties  of  methanol  at  25  °C  used  for  simulations 
are  provided  in  Table  2. 

5.  Results  and  discussion 

The  velocity  vectors  in  the  straight  portion  of  the  microchan- 
nel,  at  a  distance  of  150  p>m  from  the  channel  inlet,  is  shown 
in  Fig.  2.  This  represents  the  flow  conditions  under  constant 
zeta  potential  applied  at  both  walls  (Case  1).  The  color  indicates 
the  velocity  magnitude  +  v2  +  w2)  while  the  vectors  in- 


Fig.  2.  Flow  pattern  in  the  channel  cross  section  at  150  |jim  from  channel  inlet  (section  AA  in  Fig.  1)  for  Case  1. 
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VelocityMagnitude  -  ms'1 
0.001496 


0.0014 

0.0012 

0.001 
0.0008  H 
0.0006 

0.0004 

0.0002 

0 


1 


Fig.  3.  Flow  pattern  in  the  channel  cross  section  at  150  |xm  from  channel  inlet  (section  AA  in  Fig.  1)  for  Case  2. 


dicate  direction  of  in-plane  component  (\Zu2  +  w2).  It  can  be 
seen  that  due  to  electroosmotic  driving  force,  a  secondary  vortex 
pair  has  been  generated  in  the  flow  cross  section  such  that  the 
fluid  near  the  side  walls  moves  towards  the  bottom  plate.  The 
fluid  again  rises  in  the  central  portion  of  the  channel  to  main¬ 
tain  the  continuity.  Thus  two  distinct  vortices  are  formed  due  to 
electroosmosis.  These  kind  of  vortices  are  absent  in  case  of  flow 
without  electroosmosis.  This  modified  flow  profile  gives  rise 
to  additional  pressure  drop  in  case  of  flow  with  electroosmosis 
[12]. 

Figs.  3  and  4  show  the  velocity  patterns  for  Cases  2  and  3, 
respectively.  It  is  observed  that  due  to  the  stepwise  zeta  potential 
profile  applied  on  the  right  wall,  the  center  of  the  right  vortex  is 
shifted  upwards  in  Fig.  3  compared  to  Fig.  2.  This  is  because  for 


Case  2,  the  maximum  magnitude  of  zeta  potential  occurs  in  the 
top  region  (—90  mV)  while  in  the  bottom  region  it  is  minimum 
(—10  mV).  Hence,  the  top  region  dominates  the  electroosmosis- 
induced  secondary  flow  and  the  vortex  center  is  shifted  towards 
the  top.  Similar  behavior  is  observed  in  Case  3  where  stepwise 
zeta  potential  is  applied  on  the  left  wall. 

It  is  interesting  to  observe  the  flow  pattern  for  Case  4, 
which  is  shown  in  Fig.  5,  where  stepwise  zeta  potential  is  ap¬ 
plied  on  both  the  side  walls.  Under  such  condition,  the  cen¬ 
ters  of  both  the  secondary  vortices  are  shifted  towards  the  top 
wall. 

When  the  zeta  potential  is  linearly  varied  for  each  wall,  as 
shown  in  Figs.  6  and  7,  there  is  a  gradual  change  in  the  elec¬ 
troosmotic  driven  force  on  the  fluid.  This  results  in  shifting  of 


Fig.  4.  Flow  pattern  in  the  channel  cross  section  at  150  |xm  from  channel  inlet  (section  AA  in  Fig.  1)  for  Case  3. 


160 


A.A.  Saha  et  al.  /  Journal  of  Power  Sources  164  (2007)  154-165 


Fig.  6.  Flow  pattern  in  the  channel  cross  section  at  150  jxm  from  channel  inlet  (section  AA  in  Fig.  1)  for  Case  5. 
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Fig.  7.  Flow  pattern  in  the  channel  cross  section  at  150  fxm  from  channel  inlet  (section  AA  in  Fig.  1)  for  Case  6. 
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Fig.  8.  Flow  pattern  in  the  channel  cross  section  at  150  |xm  from  channel  inlet  (section  AA  in  Fig.  1)  for  Case  7. 
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Fig.  9.  Streamline  plots  in  the  channel  cross  section  at  150  |xm  from  channel  inlet  (section  AA  in  Fig.  1)  for  Cases  1-7:  (a)  Case  1;  (b)  Case  2;  (c)  Case  3;  (d)  Case 
4;  (e)  Case  5;  (f)  Case  6;  (g)  Case  7. 
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Fig.  9.  ( Continued ). 


the  secondary  vortex  by  a  smaller  distance  as  compared  to  Cases 
2-4.  Fig.  7  shows  the  velocity  profile  for  Case  6,  in  which  zeta 
potential  is  linearly  varied  along  the  left  wall  and  in  Fig.  8,  the 
zeta  potential  is  linearly  varied  along  both  the  side  walls. 

The  streamline  plots  for  all  Cases  1-7  are  shown  in  Fig.  9(a- 
g).  The  streamline  plots  depict  the  formation  of  vortices  for  all 
the  cases,  and  the  effect  of  variable  zeta  potential  is  clearly  seen 
in  these  figures. 

Fig.  10  shows  the  variation  of  v  component  of  velocity  along 
the  channel  width  at  z  =  50  p,m  for  the  cross  section  AA  under 
various  zeta  potential  conditions.  It  is  observed  here  that  for  Case 
1,  the  u  component  is  zero  as  the  secondary  vortices  appear 
at  the  center  line  and  results  in  a  zero  velocity  at  the  central 
core.  For  Cases  2  and  3,  the  u  component  reaches  maximum 
near  the  vortex  core  and  then  it  drops  to  zero  towards  the  outer 
edge  of  the  vortex.  It  is  to  be  noted  that  due  to  the  shifting  of 
the  vortex  core,  the  location  of  the  maxima  is  at  two  different 
positions  along  the  channel  width.  However,  when  both  walls 


are  subjected  to  the  stepwise  variable  zeta  potential,  as  in  Case 
4,  the  u  component  has  two  maxima  near  each  vortex  core,  and 
the  profile  is  bounded  between  the  velocity  profiles  for  Case  2 
and  3,  respectively.  Similar  behavior  is  observed  for  Cases  5 
and  6,  where  linear  variation  of  zeta  potential  is  applied  at  the 
each  wall  and  for  Case  7,  where  both  walls  are  subjected  to  a 
variable  zeta  potential.  It  is  to  be  noted  that  linear  variation  of 
zeta  potential  results  in  lower  electroosmotic  velocity,  which  is 
observed  in  Fig.  10,  where  the  magnitude  of  u  component  is 
lower  for  Cases  5-7  compared  to  Cases  2-4. 

The  y  component  of  velocity  is  plotted  in  Fig.  1 1  along  the 
channel  width  at  z  =  50  jim  for  the  cross  section  AA.  It  is  ob¬ 
served  that  the  velocity  profile  corresponds  to  a  fully  developed 
profile  of  steady  laminar  flows  in  square  ducts  for  all  the  cases. 
The  y  component  is  predominantly  due  to  the  pressure  driven 
flow  through  the  microchannel  and  the  electroosmosis  contribu¬ 
tion  to  this  velocity  component  is  quite  small  and  not  discernible 
in  the  figure. 
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Fig.  10.  Variation  of  x  component  of  velocity  at  z  =  50  |xm  along  the  channel 
width  at  section  AA. 

Fig.  12  shows  the  z  component  of  velocity  along  the  channel 
width  at  z  =  50  pan.  It  is  observed  that  the  z  component  reaches 
a  maximum  at  a  small  distance  from  the  wall  where  the  elec- 
troosmotic  velocity  is  maximum.  It  then  drops  to  zero  and  again 
reaches  a  maximum  value  in  opposite  direction  at  the  center  of 
the  channel.  The  variation  of  z  component  is  similar  for  all  the 
cases.  This  is  in  parallel  with  the  velocity  vector  observed  ear¬ 
lier,  where  the  flow  reverses  and  reaches  its  maximum  value  at 
the  center  of  the  channel. 

The  flow  behavior  at  the  semi-circular  bend  for  section  BB  in 
Fig.  1  has  also  been  studied  for  four  different  cases.  Fig.  13(a-d) 
shows  the  streamline  plots  at  the  bend  for  pressure  driven  flow 
and  the  applied  zeta  potential  of  —0.1,  —1  and  —50  mV  (i.e., 
Case  1),  respectively.  Fig.  13(a)  shows  the  formation  of  vor¬ 
tices  purely  due  to  the  pressure-driven  flow  with  no  contribution 


Fig.  11.  Variation  of  y  component  of  velocity  at  z  =  50  |xm  along  the  channel 
width  at  section  AA. 


Fig.  12.  Variation  of  z  component  of  velocity  at  z  =  50  |xm  along  the  channel 
width  at  section  AA. 


from  the  electroosmotic  effect.  As  the  applied  zeta  potential  is 
increased  at  the  semi-circular  bend  (section  BB  in  Fig.  1),  the 
electroosmotic  forces  cause  the  vortices  to  shift  in  a  clockwise 
direction  as  shown  in  Fig.  13(d).  The  observed  shift  of  secondary 
vortices  in  a  clockwise  direction  with  an  increase  in  the  applied 
zeta  potential  presents  an  excellent  opportunity  for  obtaining 
flow  control  in  the  serpentine  microchannel  by  regulating  the 
zeta  potential  at  the  bend.  Formation  of  additional  vortices  at 
the  top  left  and  bottom  right  comers  has  also  been  observed 
when  the  applied  zeta  potential  is  —0.1  mV,  as  shown  in  Fig. 
13(b). 

Fig.  14  shows  the  z  component  of  velocity  along  the  channel 
width  at  z  =  50  p,m  for  the  cross  section  BB.  It  is  observed  here 
that  for  a  pressure-driven  flow  only,  the  z  component  velocity 
is  zero  as  the  secondary  vortices  appear  at  the  center  line  and 
results  in  a  zero  velocity  at  the  central  core.  It  is  found  that  the  z 
component  reaches  a  maximum  at  a  small  distance  from  the  wall 
when  zeta  potential  of  —0.1,  —1  and  —50  mV  is  applied  on  the 
side  walls.  The  maximum  velocity  is  reached  for  the  —50  mV  of 
applied  zeta  potential  and  decreases  in  magnitude  as  the  applied 
zeta  potential  is  reduced.  The  velocity  component  then  drops  to 
zero  and  again  reaches  a  maximum  value  in  opposite  direction 
at  the  center  of  the  channel.  Hence,  it  can  be  concluded  that 
the  z  component  velocity  influences  the  formation  of  vortices  in 
the  semi-circular  bend  at  the  cross  section  BB  and  results  in  the 
shifting  of  the  vortices  in  the  clockwise  direction. 

It  is  to  be  noted  that  the  x  and  y  components  of  velocity  profile 
correspond  to  a  fully  developed  velocity  profile  of  steady  laminar 
flows  in  square  ducts  for  all  the  cases  in  the  semi-circular  bend. 
The  velocity  components  are  predominantly  due  to  the  pressure- 
driven  flow  through  the  microchannel  and  the  electroosmotic 
contribution  is  quite  negligible.  However,  an  offset  is  observed 
in  the  peak  velocity  from  the  central  axis  towards  the  inner  wall, 
which  is  due  to  the  centrifugal  forces  dominating  at  the  cross 
section  of  the  bend. 
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Fig.  13.  Streamline  plots  in  the  channel  cross  section  at  the  bend  (section  BB  in  Fig.  1)  for  different  zeta  potential  values:  (a)  pressure  driven  flow;  (b)  0.1  mV;  (c) 
1  mV;  (d)  50  mV. 


Fig.  14.  Variation  of  z  component  of  velocity  at  z  =  50  fxm  along  the  channel 
width  at  section  BB. 


6.  Conclusions 

A  three-dimensional  numerical  simulation  of  flows  through 
serpentine  microchannels  has  been  carried  out  for  the  channel 
side  walls  subjected  to  variable  zeta  potentials.  Different 
boundary  conditions  with  respect  to  the  zeta  potential  -  viz., 
stepwise  and  linear  profiles,  are  applied  at  the  side  walls.  It  is 
observed  that  vortices  are  developed  at  the  straight  portion  of  the 
microchannel  due  to  the  electroosmosis  and  the  characteristics 
of  the  vortices  are  quite  different  as  compared  to  a  constant  zeta 
potential  condition.  The  core  of  the  vortices  shifts  when  variable 
zeta  potentials  are  applied  at  both  the  side  walls  of  the  channel.  It 
is  also  observed  that  the  magnitude  of  the  electroosmosis-driven 
flow  becomes  smaller  for  the  linear  variation  of  zeta  potential. 
Flow  control  in  the  serpentine  microchannel  by  regulating 
the  zeta  potential  at  the  bend  has  also  been  demonstrated. 
Formation  of  additional  vortices  at  the  top  left  and  bottom 
right  corners  of  the  serpentine  microchannel  has  been  observed 
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at  the  semi-circular  bend  for  very  low  values  of  applied  zeta 
potential. 
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